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The emitted energy and the cross-section of classical scalar bremsstrahlung in massive particle 
collisions in D = 4 + d— dimensional Minkowski space Md a-s well as in the brane world M4 x T'' 
' is computed to leading ultra-relativistic order. The particles are taken to interact in the first case 

^Nj , via the exchange of a bulk massless scalar field $ and in the second with an additional massless 

scalar (p confined together with the particles on the brane. Energy is emitted as $ radiation in 
^ ' the bulk and/or (j) radiation on the brane. In contrast to the quantum Born approximation, the 

Q classical result is unambiguous and valid in a kinematical region which is also specified. For D — 4 
the results are in agreement with corresponding expressions in classical electrodynamics. 
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I. INTRODUCTION AND RESULTS 

There is little doubt that new physics will soon be discovered at the TeV scale. Absence of a light fundamental 
higgs boson implies new physics and, furthermore, a light higgs boson is unstable without new physics, of which 
examples are supersymmetry broken at the TeV scale, or maybe some variant of gauge interaction becoming 
strong at the TeV scale [ll]. 

A less well founded, but arguably more fascinating scenario for LHC would be the discovery of TeV scale 
gravity and Large Extra Dimensions (LEDs) [3| • Although extra dimensions were proposed a century ago, their 
raison d' etre, as well as their number, size and topology vary in different proposals. The potentially most 
natural justification of the existence of extra dimensions with a definite proposal for their number is string 
theory, while their most natural size is thought to be the Planck length 0(10"^^ TeV~^). Nevertheless, the 
idea that some extra dimensions could be of order 0{TeV~^) appeared in the context of heterotic string theory, 
when the size of extra dimensions was connected to the supersymmetry breakingscale Q. Earlier thoughts of 
the Universe as a topological defect embedded in a higher-dimensional space-time [J] and the development of the 
Brane- World idea led to a variety of alternative scenaria [1, Q reviewed in Q • They differ in the topology, size 
and geometry of the extra dimensions, as well as in the nature of the degrees of freedom which may propagate 
in the bulk. The technically simplest such model is the ADD 0, Q , with the Standard Model particles confined 
on the brane and only gravity allowed to propagate in the whole A^4 x T'^ space-time. 

The invention of ways to detect and measure the size of extra dimensions has attracted considerable attention 
and the study of radiation is an obvious candidate. Indeed, it has been argued that in the ADD model graviton 
bremsstrahlung from neutron collisions may give significant contribution to the cooling of supernovae [lo| . 
thus providing strong bounds on relevant parameters. Also, production of Kaluza-Klein (KK) gravitons in 
collisions of relativistic particles at LEP or LHC was extensively discussed in the literature, with emphasis 
on processes efficient to at least demonstrate their very existence [ill, [13 ■ However, the actual quantum 
perturbative calculations in the context of ADD are ambiguous, as they are plugged with tree-level divergences, 
associated with the emission of infinite momentum of modes into the compact dimensions Q . This problem is 
manifest already in the calculation of any elastic scattering amplitude in such spaces and was treated in rather 
ad-hoc ways. Contrary to the above situation, a purely classical computation, which by itself carries a natural 
mechanism to cut-off the energetic KK excitations, leads to unambiguous results [l^. Theoretically, such an 
approach is justified by the fact that at transplanckian energies scattering processes are dominated by gravity, 
which in addition, may be treated classically, at least in some range of momentum transfer [T^ . Furthermore, 
the classical formula for the elastic scattering in ADD was shown to coincide with the one obtained in the non- 
perturbative eikonal approximation, given by the summation of an infinite number of ladder graphs in quantum 
theory [13]. 

Encouraged by the above result, one may try to deal in the same way with the phenomenon of ultra-relativistic 
gravitational bremsstrahlung, and the present paper is a first step in this direction. To further motivate this 
approach, it is well known that in 4-dimensions the classical treatment of relativistic bremsstrahlung corresponds 
to non-perturbative calculations in quantum electrodynamics. Specifically, bremsstrahlung in ultra-relativistic 
scattering off a Coulomb center has been computed using the exact Coulomb solutions of the Dirac equation (see 
[l5| and references therein) . In the ultra-relativistic limit and with momentum transfer not exceeding the particle 
mass, a good approximation to the corresponding cross-section is given by the simplified Furry-Sommerfeld- 
Maue formula, which for photon energy much smaller than E"^ /m, where E and m are the particle energy and 
mass respectively, reduces to the classical result [31 ■ The result is clearly non-perturbative in the normal sense 
of expansion in powers of the fine structure constant. No analogous non-perturbative quantum calculations are 
available for bremsstrahlung in higher dimensions, especially in the presence of compact LEDs. Although certain 
aspects of classical radiation in I?— dimensional Minkowski space were discussed recently [13, 13) radiation in 
the presence of compact LEDs has not been considered so far. So, the hereby proposed classical treatment, 
which as explained in [13.] is unambiguous and expected to be trustable in the ultra-relativistic limit and small 
scattering angle, is also, to the best of our knowledge, the first actual computation of the phenomenon in higher 
dimensional space-times with the topology of 7W4 x T'^. 

The purpose of this paper is to present a classical computation of bremsstrahlung in particle collisions on 
a flat 3-brane embedded in flat space-time of arbitrary dimension and with the transverse space taken either 
euclidean (relevant to the scenario of Universal Extra Dimensions [^) or toroidal T"* (as in ADD) with equal 
radii. As a first step, the bulk interaction will be modeled by the exchange of a massless scalar field while an 
additional massless scalar </>, confined on the brane, will be introduced to mimic, in the case of toroidal extra 
dimensions, the Standard Model interactions between the colliding particles. The case of gravity differs, of 
course, crucially from massless scalar exchange. Nevertheless, the present analysis captures several important 
technical points of gravitational bremsstrahlung jT9l| , to be presented in full detail elsewhere [20j . 

Straightforward classical perturbation theory, developed in Section II, is applied to iteratively solve the 
particle equations of motion and the field equations as well. The radiation efficiency e, i.e. the fraction of the 
initial energy emitted as bremsstrahlung radiation, is computed in Section II to leading ultra-relativistic order 
and low momentum transfer and is found to be 

s,^%L^C, (7^3/63)1+^, (I.l) 
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where is the Z?— dimensional classical radius of the colliding particles and b the impact parameter. The 
constant Cd is given in (|TL47|. 

The case of toroidal extra dimensions is dealt with in Section III. The light <i>— KK modes participate crucially 
both in the interaction between the colliding particles and in the radiation emitted in the bulk. These roles 
are studied separately with appropriate choices of their couplings to the radiating particles. The radiation 
efficiency, depending on which interaction is assumed to be dominant (via $ or 0) and on the nature of radiation 
(i>— radiation in the bulk or (p on the brane), is to leading ultra- relativistic order shown in Table ID Note 
the interrelations among the various entries of the Table as well as with (jl.ip . The range of validity of the 
computation, the consequences of the results and the modifications expected in actual gravity are presented in 
the discussion Section IV [13]. Finally, three Appendices in the end contain derivations of formulas used and 
clarifications of technical approximations made in the text. They also contain a qualitative explanation of the 
basic features of Table I. 



Radiation 

Interaction 


Bulk M4 X T"^ 


Brane Ait 


Bulk M4 X T'' 






Brane M4 







TABLE I: Radiation emission efficiency in M4 x T''. The standard "electromagnetic" (-Ro) and the D— dimensional 
"gravitational" (r^) classical radii of the colliding particles, as well as the constants Cab and Co are given in the text. 



II. SCALAR ULTRA-RELATIVISTIC BREMSSTRAHLUNG IN Md 



In this section we formulate our approach in the simplest case of two massive ^ scalar point charges interacting 
via a linear scalar field in Minkowski space of arbitrary dimension D. This will serve as a reference calculation 
for the subsequent treatment of compactified extra dimensions. 



A. The model - Generalities 



Consider two particles with masses m and m' moving along the world- lines = z'^{t), x'^ = z"^{t') 
interacting with a massless scalar field $ with coupling constants / and /', respectively. The relevant action is 

S = - j dTyf^[m + f^{z)]~ j dr'\/i^[m' + /'$(z')] + ^ j d^xdn^d'^^^, (II.I) 

with M, A, B, . . . = 0,1, . . . , D — 1 — 3 + d, and the metric signature is (+, —,—,...,—). P = ijmnz^^ z'^ , and 
the dot denotes differentiation with respect to r or r'. The classical radius of the particle m is defined by 
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m 



— ' d ' 



(II.2) 



and similarly for to'. 

Linearity of the field equations implies that $ = $,„ + , where 



□ = p{x), □ = p'{x), □ = T]^^^dMdN, (II.3) 

the sources being 

p(x) = ff dT^S^ix - z{t)), p'{x) =ff dr'^d^ix - z'(t')). (II.4) 



The interest here is in strictly massive particles. Proper treatment of the massless particle case would require the use of the 
Polyakov action for them. 
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The particles' equations of motion read 

^ [(m + = (II.5) 

and similarly for m! , provided the parameters t, t' are chosen so that 

i(r)2 = 1 = i'(r')'. (II.6) 

This can be rewritten as 

(m + = / n^^^ajv$, n*^^ = r^*^^ - i^^i^, (lu) 

where H''^^ is the projector onto the space orthogonal to the world-line. In this equation the field on the right 
hand side contains both the self-action term and the mutual interaction term The first gives rise to 

divergences, taken care of by mass renormalization, and also for D ^ A to the introduction of higher derivative 
classical counter-terms JT^, which for simplicity will be ignored. Thus, the equations of motion become: 

Fourier transform using ^'(a;) = / '^{k)e^^^'^d^k/{2TT)^ and '^{k) = J '^{x)e^^'^d^x. The retarded solutions 
of the wave equations pi. 31) can then be written algebraically: 

$(fc) = G,et(fc)p(fc), G,et(fc) -G(fc) + z^£(fc")5(fc2), Gik)^-V^, (11.9) 

with the Fourier-transforms of the source terms given by: 

p{k) ^ f j e'^-'-^^Ur, p'{k) ^ f j e^'^-^'^^'^dr'. (11.10) 

The energy- momentum P^^ of the scalar radiation, emitted by the particles during the scattering process, 
can be computed by integrating the field stress-tensor 

between the space-like hypersurfaces S±oo corresponding to i — > ±oo and reads: 

P^' = j T'^^dSN - J T^'^'dSN. (11.12) 

Given the fall-off of the stress-tensor at spatial infinity T*^^ ^ l/r^^^-f ©(r^*^^^^)) (derived from the behavior 
of the Lienard-Wiechert potentials), one can write P*^ as an integral over the closed boundary of the space-time 
tube and, subsequently, transform it to the volume integral P^^ ~ J d^xW nT^^^ , which upon substitution of 
(|IlTT1) leads to 

P^' = J d''x{p + p')d^<!>, {11.13) 

or equivalently 

^^ = -(^ j d^kk^'m[p*{k) + p'*{k)]. (11.14) 

Substituting here the retarded solution (III.9[) and taking into account that the contribution of the principal part 
G{k) — ~ V (fc^)^^ vanishes being odd under fc*^ — > — fc^^, one is led to 

= j d''ke{k^)k'H{k^)\p{k)+p\k)\\ (11.15) 

Finally, denote the frequency = oj, use 6{k^)5{k'^) = 5{u) — |k|)/2|k|, and integrate over |k| to obtain the 
spectral-angular distribution of radiation: 

liP'^^ 1 



dujdn 2(27r)^- 



^L.^-'V^p{k)+p'{k)\' (ILlf 



fe*f=Lj(l,n) 



where or n parametrize the unit sphere in the {D — 1)— dimensional Euclidean subspace. 
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B. Perturbative solution 

1. Particle trajectories 

We next solve the equations of motion for the particles and the field $ iteratively in powers of the coupling 
constants / and /'. For the particle trajectories we write 

z^\t) = h^' + u^'t + (5z*^(t), z'^^(t) = u'^'t + 5z'^^(t), (11.17) 

(assuming z^^(O) — 6*^, z'*^(0) = 0) where u^'^ and u'*^ are the unperturbed constant four- velocities, specified by 
the initial condition and chosen to satisfy = 1 = u'^ . Combine with (jll.6|) to conclude that the perturbations 
of velocities must be orthogonal to the unperturbed world lines: 



UA^i6z'') = 0, u'A^iSz'^) = . (11.18) 

To zeroth order, the solution of the wave equation describes the non-radiative Coulomb field of two non- 
interacting particles, which in terms of the Fourier-transform reads 

<i>m{q)^-^e''-'^°^Siq-u). (11.19) 

Substituting this into the equations of motion (|IL8p we find for the first order perturbation: 

^^(t) = ^n^^dB'fm'ix) ^ ^ , (11.20) 



dr2 



where 11"^^ — rj^^ — u^u^ is the projector on its unperturbed world-line. The solution for 5z is 

= -» 1 / (^'^^'"^^ -l^iq-ur) 5{q ■ u')U''''q^ , (11.21) 

with an analogous expression for 5z' and with the exchange w u', & — > —6. They automatically satisfy the 
orthogonality conditions (|II.18|) . 

2. The radiation and its main features 

The next step is to compute corrections to the field $ due to these perturbations of the trajectories. These 
will correspond to the lowest order radiation field. To compute radiation using pi.l6p one needs the Fourier 
transform of the corresponding sources. Our treatment is symmetric in m and m', so we write only unprimed 
quantities. The zeroth order terms do not contribute to radiation, since the corresponding Fourier transforms 
pi.igp vanish on the light cone fc^ = (except for the trivial point =0). Ditto for the last two terms inside 
the parenthesis in the integrand of (|II.21|) . Thus, to lowest non-trivial order ^: 



p{k) c^if dre*^(°)+'('=-")^(fc • dz{T)). (11.22) 



Substitution of (III.2ip gives 



(27r)^-2 m{k-uy 



with 



The corresponding integral I'^ik) for the partner particle is obtained from I a by the substitution u ^ u' , za{0) o 
z'a{0), bA — ^ —bA- The integrals are evaluated in Appendix I in a fully covariant form. 



^ A word of caution about notation: To avoid the introduction of too many symbols, we use the same symbol p (a) for the generic 
source in the previous section, (b) for its zeroth-order and first-order expressions, (c) for the sum of the two that one has to 
substitute in the energy loss formula l lII.lSII . with the zeroth order term giving zero contribution. Hopefully, the context makes 
this notational simplification harmless. 



6 



It is convenient to work in the lab frame, in which u'^ = (1, 0, 0, 0), i.e. with m' initially at rest. Define 
the remaining kinematical variables in this frame by 



u'* = 7(l,v), ~ LL){l,n), k-u — ^u}^, ^=\—vcos6, k-u' — uj, 



(11.25) 



where as usual 7 ~ ~ v'^ and 6 is the angle between n and v. The invariant expression for the Lorentz 

factor is 7 = u • u'. 

Furthermore, we define the zeroes of the parameters t and r' in such a way that at the point of closest 
proximity of the scattered particles they satisfy u ■ b — ~ u' ■ b. In other words, the "initial conditions" 
defined at r = = r' are given at the moment of least distance between the particles. But then, b'^ — and 
V • b = (easiest visualized in the CM frame). The vectors v and b define the collision plane, while the generic 
{D — l)-dimensional unit vector n can be decomposed as: 



V b 

n = - — - cos 9 + - — - sin 9 sin t/j + m sin 9 cos ■(/;, 

V b 



where m is a Z? — 1 dimensional unit vector orthogonal to the collision plane. Then 

k • b = — fc • 6 = buj sin9 sinip, 
where the impact parameter b = | A | is also given by the invariant expression: 



,2 [(u'•6)u-(u•6)^i'l2^l/2 



(11.26) 



(11.27) 



(11.28) 



Occasionally, it will be convenient to express the angular dependence of radiation in the rest frame of the 
moving particle m. In that frame the radiation angle will be denoted by 9q, which is related to 9 via 



1 / 1 



cos 6*0 - , . „ 
V Vs7 

Using formula (jV.12[) of Appendix I one obtains 

j2 j/gifc-2(0) 



1 



sin 9q 



sai9 

IT' 



p{k) 



(27r)"+im73w362n^^ 

where n = D/2 — 2 = d/2. The "hatted" functions Kn{z) = z'"-Kn{z). Their invariant argument is 

(k ■ u)b buj£^ 

V 

Similarly, the leading contribution of m' to the radiation field in the lab frame is 



p'(fc) 



with 



(27r)"+im'72w^&^" z' 

, (fc • u')b buj 



J cos 9z' Kn{z') + i sin9 sinipKn+i{z') 



Substitution in pi.l6p gives the angular and frequency distribution of the emitted radiation. 



(11.29) 



(11.30) 



(11.31) 



(11.32) 



(11.33) 



3. Frequency and angular distribution. 

(a) The low-frequency limit. Using lim Kn{z) — 2"~^r(n) for n > and lim Ko{z) ~ ln(2/z) — Ce (where 

Ce is the Euler constant) in (|II.30I) one concludes that for w — > the leading term for any D will be the second 
term in (jll.30p and 

^ »/VT(n+l) sing sin 1 
' u^o~ (27r)«+im7362n+i (1 -i; cos 61)2 ■ ^ ' ' 

The divergence in w, is reminiscent of the infrared divergence of the corresponding Feynman diagrams. However, 
upon multiplication by uj^~^ from the measure and integration, it contributes a finite amount to the radiation 
loss. Furthermore, the non-zero value of dE-^g^^/ duj{uj = 0) for I? = 4 is tiny for large 7. 
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FIG. 1: The contribution to the frequency distribution of the emitted energy of the real (sohd hne) and the imaginary 
(dashed hne) part of p for D = 6 and 7 — 20. The two curves are rescaled by the same factor. 



^max — ^cr 



?=o b 



(b) Frequency cut-off and beaming. The exponential decay of p{k) due to lim Kn{z) ^ \ 772" ^^"^e ^ gives 

>oo y 2 

a 0— dependent upper cut-off for the frequency of the emitted radiation: 

o(l — vcoau) 

In the ultrarelativistic case 

- ,(^.^\-.) . (n.36) 

so that most of the radiation is beamed inside the cone 

e < 7^1 (11.37) 
in any dimension D. The maximum frequency of radiation is 

27' 



(11.38) 



On the other hand, p'{k) does not exhibit sharp anisotropy in 9 and the associated frequency cut-off in the 
ultra- relativistic case is 27 times smaller than Wmax: 



max 



1 ^ (11.39) 

27 



(c) The imaginary part of p{k) is negligible. One might expect, that as in D = 4 electrodynamics one can 
combine the low-frequency amplitude (|II.34p with this cut-off to estimate the total radiation loss. This turns 
out to be incorrect here. The imaginary part of p in pi.30p is suppressed compared to its real part by a factor 
sin 9 ~ < 1/7. So, the leading ultrarelativistic contribution to the radiation loss is due to the real part of p, 
i.e. the first term of the amplitude pi.30|) . This is demonstrated for D = 6 and 7 = 20 in Figured) 

(d) p'{k) is also negligible. Using (IVI.SP and (jVI.17p one may convince oneself that in the lab frame the 
contribution to the total radiation (a) of the particle m' and (b) of the cross-term in (|II.16p are both negligible 
in the ultra-relativistic limit. In particular. Figure [5] shows separately the contributions of m (solid line) and 
rn' (dashed line) to the emitted energy in the lab frame. So, the radiation from m' is important at very low 
frequencies, but may be neglected everywhere else, as well as in the total energy loss. Thus, we shall neglect 
the contribution of everything else but the real part of p{k). This gives most of the energy lost in the collision. 



(e) Finally, it is convenient to introduce the angle 9 by 

cos 9 = V, sm9 — 



1/7. 



(11.40) 



8 




y = 0}b/y 



FIG. 2: The relative contribution to the frequency distribution of the emitted energy of p{k) (soHd hne) and p'(fc) (dashed 
line) for D = 6 and 7 = 20. The two curves are rescaled by the same factor. 



Within the cone ^ 6 ^ 6 the quantity ^ for ultrarelativistic velocities remains small varying in the region 
1/(27^) ^ ^ ^ 1/7^- Note that for 9 = 9 the m-rest-frame radiation angle is 6*0 = tt/2. 

(f) The angular distribution dE-cad/dO- As argued above, the leading contribution to the total energy loss is 
due to the real part of the amplitude p{k), i.e. 



p{k) ~ Re p{k) 



(27r)"+imw37352n — i;C0S 



1 



(11.41) 



which does not depend on angles other than the main polar angle 9. 

Substitute (|lL4T|) into (|IlT61) and use (fVLSl) to integrate over the frequencies. The result is the angular 
distribution of radiated energy: 



dE, 



rad 



0Fr 



dn 8(27r)2i5-3r(£) - 2) 63D-9^6,„2 

Integration over gives to leading ultra-relativistic order the 0— distribution 



r^''+'^ sin^-3g. 



d9 ~ 7r3(-D/2-i)r2 

In the special case f = f and m = m' one obtains for the radiation efficiency Ed = i?iad / E 
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d9 28+3'i7r3+3''/2r2(l + d/2) 
plotted for various dimensions, 7 = 500 and rd/b = 0.3 in Figure |31 



(11.42) 



(11.43) 



(11.44) 



(g) The frequency distribution dE-ca,^/ duj . Similarly, define y = 0^6/72 and substitute (|II.4ip into (jlLlGj) to 
obtain 



d£d 



l+d/2 



__^3+2. {'^y'^'^f^'^J d9 [1 - (7^0-^^.^/2(2/7^0 Si 



sini+'^e. 



dy (27r)5+2'ir 

evaluated numerically and shown for various dimensions, 7 = 500 and rd/b = 0.3 in Figured 



(11.45) 



4- The total radiated energy Emd 



Finally, using (|VL5p one may integrate pL43p over 9 to obtain the leading ultra-relativistic limit of the 
radiation efficiency (for m = m! and f = f) 



E, 



ad 



E 



3\ 



(11.46) 
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FIG. 3: The angular distribution dsd/dO versus 6/6 for 7 = 500, rd/b — 0.3 and for d = (solid line) and d — 2 (dashed). 
The maxima are at 6(1,2 ^ ^''^ 6. 



de_, 

a 




FIG. 4: Frequency distribution dsd/dy versus y = hu/'y^ for 7 = 500, rd/b = 0.3 and for d = (solid curve) and d = 2 
(dashed line). 



with 



_ (2 + d)r(^)r(2^] 



27+2d7r3(2+d)/2p(-4 ^ ^) 



(11.47) 



The radiation loss grows with d as 7 + . For d = 



£0 



3 X 



29 7r2 ^ V 6 



in agreement (modulo polarizations) with the corresponding formula in classical electrodynamics 



5. The bremsstrahlung cross-section an 



The energy differential cross-section with dimension (length) 2+'^, i.e. the fraction of energy d{E^i^d/E) emitted 
per unit time and per unit incoming beam flux J = dN / dA-i-^-ddt is given by: 

^i+d^h , / p^ ^h+A |p(fc) |2fei+dd6 , 



dad = ni+db'+''dhd{E,^i/E) = 



2(27r)3+'^m7 
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with the volume of the unit (n — 1)— dimensional sphere being fln~i = 27r"/^/r(n/2). Its integral over b from 
^min to infinity is the frequency and angular distribution of the emitted energy fraction per unit time and unit 
incoming beam fiux, while its total integral is the fraction of initial energy emitted per unit time and unit flux. 
Using (|II.46|) this is 

/ \ 14-2(i 



7717 J l + 2d \bmm ' 



6. The average number (N) of emitted quanta 

A useful, though not classical, quantity to estimate is the number of emitted quanta. Divide the right hand 
side of (jII.16[) by the quantum energy huj to obtain 



fc"=w(l,n) 



substitute pi. 411) and integrate, as above, over frequencies and angles to obtain an estimate of the total number 
of emitted $— quanta, not counting the ones from m' with very low — (see Figure [2|) 

(^.) , rmr^(i + .)r^(^) ^ .^y.a. 

^ ' (3 + d) 25-'i7r4+3d/2r(2 + 2d)r(2 + d) h ^ \ b J ' ^ ' ' 

III. SCALAR ULTRA- RELATIVISTIC BREMSSTRAHLUNG IN X4 x T'' 
A. The model - General formalism 



The above analysis is easily extended to study bremsstrahlung in the context of a scalar multidimensional 
model with the extra dimensions compactified on a torus T''. The model of interest describes two particles on 
the brane with masses m and m', respectively, and two scalar fields. One, the "scalar graviton" $, lives in the 
D ^ A + d dimensional bulk and is supposed to imitate gravity, while the other, cj), lives on the 3-brane and 
mimics the Standard Model forces. The action of the model is 

particles 

with ^ = 0, 1, 2, 3; M = 0, 1, ...,£>— 1. z^(t) and z''^{t') are the trajectories of the particles with masses m 
and m', respectively. Notice that no $ — </> interaction is included, because it is not relevant in the leading order 
computation that follows. The coupling constants e, e' are dimensionless, while /, /' have length dimension d/2. 
With all radii of the torus taken for simplicity equal and denoted by L, the kinetic term of the zero mode of 
$ is {L'^/2) J d'^xdfj.^d'^^. It is brought to the standard normalization by the rescaling $ -> which 
converts the interaction term to J dr^/i?^ with the 4-dimensional coupling being 

/4 = fL-"^ . {111.2) 

The particles, apart from the radius defined in (|II.2|) . are also characterized by the "electromagnetic" classical 
radius, defined by 

Ro = e^/m , (III.3) 

and similarly for the particle m! ^ . 

This simplified model is rich enough to study all four cases presented in Table U namely, the energy emitted 
in the bulk or on the brane, when the scattering of the two particles is dominated by either the $ or the 
interactions. The cases of 0— interaction with (/)— radiation on the brane and interaction with $— radiation 
in the bulk are special cases of the previous section. The formalism necessary for the study of the crossed 
situations i.e. interaction via $(0) and 0(<I>)— radiation in the bulk, is the topic discussed next. 



^ In the quantum theory, one may trade / for a "Planck" scale M* writing / = {M*/h) "^Z^. Then, one obtains fi = {M,L/h) '^/'^ ^ 
the analog of the ADD relation between "Newton's" constant, "Planck" scale and internal space volume. 
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Retarded propagator for massive modes and particle trajectories. Denote the space-time coordinates by = 
{x'', y*}, with ^ = 0, 1, 2, 3; i = 1, . . . ,d and use bold letters to denote the rf— vectors or the extra components 
of the £)— vectors. 

The retarded Green's function of the d'Alembert equation 

OoGDix, x',y- y') = 6\x - x')<5'^(y - y') (111.4) 

is now expanded in Fourier series: 

Gd{x -x',y- y') = —1— / d^pe-'P<^-^'^ ^ p , (111.5) 

{2Tr)Wd J ^ p^ -n^+ iep° 

where n = 27rn/L, ri = \/fi^, and = L"^ is the volume of extra space. Correspondingly, the solution of 

an^x,y) ^ PDix,y) (111.6) 

with the source locaUzed on the brane ^ 

PD{x,y) = p{x)6''{y) (111.7) 

is 

= (2^ / ^'P^--'-' E ^Sh^o Pi^') d'^' (ni-8) 

Its restriction to the 3-brane can be rewritten using the four-dimensional propagator 

$(x) = $(x,y) = I Gi{x-x')p{x')d^x\ Gi{x-x')=GD{x-x',y-y') , (111.9) 

y=0 J y=y'=0 

whose Fourier transform is 

<^4(P) = E 2 -l^- - 

Vd ^ , »^ — n^ -h lep^ 

In the four-dimensional language the massive scalar KK modes act as independent massive fields universally 
interacting with scalar charges through the total field 

^{x) = J2 *n(a;). (III.ll) 

Thus, keeping four-dimensional conventions for the Fourier-transform, we can present the fields generated by 
the particle m to lowest order as 

$„(g) = -^e^«-W^(9 • u)^^, (III.12) 

and similarly for m'. Taking into account only mutual particle interactions one obtains the lowest order correc- 
tions to the particles' trajectories: 

It can be checked that the corresponding corrections in the transverse directions 
by parity and, as expected, the particles do not leave the brane. 



* In this section we have to distinguish the D-dimensional scalar source defined in the whole £> = 4 + d space and the 
four-dimensional p(x) localized on the brane. 
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^—radiation flux. Consider a {D — l)-dimensional space-like hypersurface in the Z?— dimensional space-time 
and choose it for simplicity to be orthogonal to the time axis. The total momentum of the field $(a;,y) 
associated with it reads: 

P^' = f T^'^d^'xdy. (III. 14) 

The field generated by a source localized on the brane is 



nS 



where p{k) is the four-dimensional Fourier-transform of the source. Substituting this into (jIII.14p and integrating 
over the three space and the torus, one obtains for M — i: 

Vit) ^ I dk^dk- V P(fe".k)p-(fe-o,k)e-(^"-^-Vfc"n 

J J^, [(P)2-k2-n2][(fc'0)2_k2_n2]' 

where k is a wave 3-vector on the brane. The sum vanishes by parity and reflects the expected fact that the 
radiation is emitted in the bulk symmetrically with respect to the brane. 

Similarly, the change of the tangential to the brane components of the momentum of $ between two space-like 
hypersurfaces I]±oo of topology x T'' is again given by pi.l2p . or equivalently 

= /" d'^x [ dyVjvT^^, (III.17) 

from which 

= y d^x dy {pD + p'd) 9'^* , (111.18) 

with the sources localized on the brane. Thus, the integral reduces to a four-dimensional one, and the retarded 
scalar held can be computed using the four-dimensional Green's function. In particular, the energy emitted is 

oo 

where uj = |k| and p^°^{k) is the four-dimensional Fourier-transform of the source. 



B. Particles interact via <li and emit (^—radiation on the brane 

Consider first the case in which the colliding particles interact mainly via Technically one may take e' = 
in ()III.1|) . To compute the 0— radiation emitted on the brane by particle m, one has to solve the field equation 
for (j) with the source term of order //'. The radiation amplitude is the Fourier-transform of the source p{k) 
on the mass-shell fc^ = of (p, while the radiation energy loss will be given by the n = term in (jIII.19[) . or 
equivalently, by (|II.16p with D = A and p' = 0. 



1. The emitted energy 

The source p{k) is of the form (|II.23p . with an obvious change in the couplings in front and summed over all 
KK modes. Using (jV.18IV.19p relevant to the contribution of the generic massive mode n, one obtains 

p{k) = 27:mLd^'^v'^^^ E [bzKo{Z,,){^k ■u-k-u')+ ivki{Z^)k ■ &) , (III.20) 
with Zn = \/ -V v?}?' . 

The sum over the massive modes depends on the ratio hjL. For small b/L a, large number of modes will be 
excited and for a summand that depends only on the magnitude of il one can write (see Appendix III for an 
estimate of the error due to this approximation) 

$]/n « y /„dn = -^^d^i j hn'^-'dn. (III.21) 
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FIG. 5: Frequency distribution dsd/dy (d=l, solid curve) and ^de/dy (d=2, dashed) versus y = 61^/7^, both for 7 = 500, 
rd/b = Ro/b = 0.3. 



Define the new variable w =1 + {nb/ zY and use formula ^il] (for a, /i > 0) 



(III.22) 



to obtain 



(III.23) 



where Ku{y) = K^^{y). The resulting expression for the source coincides with the (4 + (i)-dimensional (jll.SOp . 
According to the reasoning of Section II, the leading contribution to the emitted energy in the ultrarelativistic 
case comes from the first term inside the parenthesis of (jIII.20p and gives 



eff 



(III.24) 



Thus, with an appropriate correspondence of the couplings, the source looks the same as for the massless field 
in the case of non-compactified space-time of dimension D — A + d. The total emitted energy, however, will be 
different because of the different phase space in the two cases. 

Substitution of (jIII.24[) into (IIII.19P and integration over the frequencies leads to the angular distribution 



dE, 



rad 



dn 25(27r)4+rf r (^) m27663-h2d 
or equivalently (for J — f and m — m!) 



dVL 



25(27r)4+'ir(^) b 



(III.25) 



Thus, up to the overall coefficient the angular distribution is the same as in AD shown in Figure [31 No sign 
of extra dimensions in the angular profile of the radiation emitted on the branc, in contrast to the frequency 
distribution shown in Figure [5] 
Finally, use 



j d^ j Sine CHl- 1^0^ d6 



An 



7^ + 0(7'), 





to integrate over the angles and obtain for the efficiency in the ultra-relativistic limit: 

Ro /rrf\2(H-<i) 



7, 



(III.26) 
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with 

r(2M)r(3^; 



(III.27) 



3-24(27r)3+rfr(^) ■ 
As a check, notice that for d = and the identification e = / it coincides with pL46p . 

2. The cross-section 

For the scattering process, which takes place on the 3-brane, the total cross-section for ultra-relativistic 
scattering is the integral of the differential energy cross-section 

rfcr$0 — 2TrbdbdEyg^(i/mj, (III. 28) 

with d£'iad(w, fl, b) given above. Integration over the energy uj of the emitted scalar and the angles fl — {9, ip) 
on the brane leads to £'rad(&) and upon integration over b to (for f = f and m — m!) 

a^, ^-^C^,rl^({±-T 7. (111.29) 
For d = it coincides with the corresponding expression derived in Section II. 

C. Interaction via (j) on the brane, emission of $ in the bulk 

Consider the action (|III.1I) with f — 0. Now the particles interact only via the brane field (j), but radiate 
both on the brane (</) and massless mode of $) and in the bulk (massive modes N of $). One wishes to estimate 
the amount of radiation into the bulk. 



1. The energy radiated in the bulk 



Following steps analogous to the ones above, one starts with the source for the A^— th mode (eqn. (jII.30P with 

D = 4) 



ee'f 



2tT L'^ 7717^ 



7z 



7 ] Kq{z) + iv sinO siml}Ki{z)/ z 



(III.30) 



with the arguments of the Macdonald functions in p{k) and p'{k) being 



k ■ ub b 



- _|_ ^2 _ (^og 



and z' 



k-u'b 



■VuPTTp, (III.31) 



depending on the mode vector N. 

Next, one has to substitute pil.30|) into pil.191) . integrate over frequencies and angles and, finally, sum over 
the KK tower. 

As in the previous cases, the main contribution to the emitted energy is due to the real part of the fast 
particle's source. To compute it, introduce S = fc • u/{'yk ■ u') and write 



E, 



ad 



2(27r)4+'im276 



NS 



(III.32) 



Isolate the contribution of the massless mode and use (jIII.2ip to convert to integration over masses. Then change 
to polar coordinates (w, a) hy uj — wcosa, N = wsiaa with < a < 7r/2 and defining S = 1 — iicosacos^, 
integrate over w: 



d+2j^2,r.^,. {v/bE)d+^2^,fd + 3 

Jo 



r(d + 3) 

Then the leading in 7 contribution to the total emission in the bulk is given by 

7r/2 

Erad = F I cos^ a sin'*^^ ada 



dO sin e 



"d+5 



1-7^2)^ , with F = 



mj^ J (27r)4+'ir((i-t-3)6'^+3' 



(III.33) 
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The integral over 9 can be converted into an integral over S from = 1 — u cos a to = 1 + w cos a. The range 
of integration over a is such that one can omit the terms with ^q, and write 

""'^ /,-(d+4) .-(d+3) .-(d+2)\ 

E„,-.F I c„,asi„-.a.„ f^-2.^«^+yf^ , ,,,1.34) 

^ ^ 

with all remaining integrals of the form (jVI.Sp of the first kind {2m > n + 1). The difference in the range 
of integration is insignificant, because it leads to 0(7"^) contribution in the integrals, which is negligible in 
the ultra- relativistic limit. Change the range of a integration to [0,7r] and use (jVI.51) to evaluate piL34l) . To 
leading order in 7 the result is (for e = e' and m = m!) 

e,^ {ee'Ur [j) ^ ^ C,^ ^ {^) (III.35) 



with 



p3 



2. The cross-section a^if. 



Upon multiplication of piL35p by 27rfe db and integration over h one obtains the total cross section 



D. Interaction via $ and $— radiation in the bulk. 



Consider, finally, the case in which the particles interact via exchange of $ and emit radiation in the bulk. 
Their couplings to $ are / and /', respectively. The source of radiation of the N— th mode is then: 



with the argument of the Macdonald functions being 



, Lo sin 6 sin ■0 
iv" , 

+ N^z 



{y uj"^ + N"^ - Lov COS 6^ j +\? 



1/2 



(III.37) 



(III.38) 



and the dimensionless products fc • ub/^v and k ■ u'b/^v given by pil.31|) . 

Assuming again a large number of interaction modes, we replace the sum over modes by integration according 
to (Hirni) : 



(27r)''/2+iTO73i,3fod 



2 uj sin 9 sin tp - 

V W + N Z 



(III.39) 



Thus, the argument of the Macdonald function corresponds to the massless interaction mode pil.30|) . so again 
the real part gives the main contribution to the emitted energy, which in polar coordinates (w, a) in the oj — N 
plane, becomes 



p(^)(fc) = - 



(S-i-7')^d/2(2), ^=l~v 



(27r)d/2+i^y„3^d 
Substitute it in equation piI.19|) and integrate over w, using: 



cos a cos ( 



(III.40) 



dww^+^kl/^iwb^) 



4r(d-|-2)(6S 



Perform the remaining integrations over the angles Lp, 9 and a and use similar approximations to obtain for the 
leading ultra-relativistic contribution to the emitted energy (/ = f , m = m!) 



~ Cd ftf'\ ( I 



3<i 



L\ 7 



3+3d 



(III.41) 
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with Cd given in pi.47l) . As expected, this expression is identical to the one obtained in the case of 
D— dimensional Minkowski space (|II.46|) . even though the intermediate formulae and numerical coefficients 
in the phase space integrals are different, corresponding to different spatial topologies. This is a consequence of 
our approximation to convert mode summation to integration. 
The corresponding energy cross-section is 

n / Y \ 3d „ / \ l+3d 



IV. DISCUSSION - VALIDITY OF THE APPROXIMATION - PROSPECTS 

An unambiguous classical computation of bremsstrahlung radiation in ultra-relativistic massive-particle col- 
lisions was presented in the context of a simplified scalar model in arbitrary dimensions, of which some may be 
compact. Scalar fields were used to model both the graviton in the bulk and the standard model interactions 
on the brane and the main results for the radiation efficiency are summarized in (|I.ip and Table I. A quick 
glance at these leads to the following remarks: (a) In all cases there is enhancement of the efficiency by factors 
of 7. (b) Radiation emitted on the 3-brane is enhanced by one power of 7, while (c) each compact large extra 
dimension to which radiation can flow, contributes to the efficiency an extra power of 7. 

A few comments are in order concerning the validity of our approximations. The leading order perturbative 
classical computation per se is a good approximation as long as the interaction energy is much smaller than 
the total available energy, i.e. ^ C/(6) ~ <^ E c± mj. On the other hand, the classical approach 

is relevant if one can justify (a) dealing with particle tra ject ories, and (b) treating radiation classically. As 
has been discussed in non-relativistic quantum mechanics 111'] and applied to the relativistic case as well [2^ , 
requirement (a) implies (al) small angle scattering and (a2) |[/'(6)|6^ ^ ft. Therefore, the impact parameter 
has to satisfy 

«6«^ =b,. (IV.l) 



In addition, condition (b) of classicality of the radiation is written as 

iN)-j-[,^) -J-ed»l. (IV.2) 

Notice that for ^ 0(1) the latter is satisfied for m6 ft7, which is equivalent to ftwmax ^ h-y'^/b <^ 7717, 
which also has to be satisfied, since an emitted quantum cannot carry more than the total available energy. 
Finally, the deflection angle in the ultrarelativistic case is given for small momentum transfers |i| ^ s by ^ 



Sa 



li- Al 



T-^ — OO 



T^oo myZiT)^ ^07 J q'^[q-u) 



|A|7 

where b = |b|. Use lim sin(g • uT)/{q ■ u) = 7r5(g • u) and perform the integration to obtain for f — f 



Sa — 



rii+d ^ b J 7 

The region of validity of this formula is 6a <C <C 1 or equivalently b ^ 7-^ 'y-3/2(i+<i)^ and follows from 

the condition of small momentum transfer stated above. 

For ft — >■ the constraints are summarized to 6 ^ rd7~^^^^'^. Taking into account also the quantum 
constraints with d 7^ the impact parameter is restricted to hj/ni <^ b <^ be, which requires (for ft = 1) 
7 ^ fl^'^mL. In the special case of d = the region of validity of our approximation is 6 3> fv) /m. 

The results of the present paper are suggestive, but by no means conclusive about gravitational bremsstrahlung 
itself. Gravity has many and significant differences from the toy scalar model, (a) In gravity there is an 
additional relevant scale Rg — Gn\/s, which in principle enters the condition for the minimum classically 
reliable value of the impact parameter, (b) In addition, the fact that gravity couples to the energy-momentum 
is expected to lead to extra enhancement of radiation in the transplanckian regime. Finally, (c) gravity is non- 
linear with important effects due to these non-linearities. These qualitative features of classical transplanckian 
gravitational bremsstrahlung have been verified and presented briefly in 19], while a longer detailed exposition 
is in preparation [20l |. 



^ The discussion here concerns the case of jMo. Similar analysis can be carried out in 7VI4 X Td- 

® Note that: (a) The cosine part of the exponential exp(—iq ■ ur) in 1III.2II 1 does not contribute in 5a. (b) As explained below, the 
points T = = t' on the trajectories are chosen so that u ■ b = 0. 
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V. APPENDIX I: MOMENTUM INTEGRALS 

(a) Start with the following invariant integral 

J J 

In the frame with u' — (1, 0, 0) one may immediately integrate over dq'^S{q'^) to obtain 



.^V(^-+q;")-"". (V.2) 

q 



Introduce 7(1, v), with 7 — 1/Vl — and decompose q and b along v and perpendicular to it, writing 

q = g||n + qj^, b = &||n + bj^, where n = u/jv = v/v. It is straightforward to check that 

bn= — {ju'-b-u-b) and |b_L| = -& - / n • (V.3) 

jv \ [u ■ u Yv"^ J 

Integrating over dq\\ with 5(k ■ u + q\\"iv) we obtain 

/ = ^-^>^-uHhv / e ^ _LJ^c-*"''ii/^'^i^„(z), (V.4) 

7W J qj^ -I- (fc • uY l^^v^ 7w6^" 

where n = — z = k-ub/'yv, Kx{w) = w^Kx{w) and b = |bj^| now. For 6|| '-^ u-b = (the case of interest 
in the main text) the above simplifies to 

^ = -^^"W- (V.5) 
Correspondingly, the primed integral obtained by a it o u' and 5m o —6m exchange, 

^ = 7 ^ 1 -2 = ^^^«(^ ) (V.6) 

with z' = k ■ u'b/^v. 

(b) Consider next the vectorial integral 



It may be computed with the help of 



.81 .( dl db\\ dl db \ 



'db^^ \db\\ db^' db db'^^ ) 
Its explicit fully covariant form is obtained using (|V.4|) and (|V.3|) 



dl_^.k^ dl ^ I Kn+ijz) 

\v ' db b K„(z) ' ^ ' 



where use was made of the formulae 



71 n 

KXz) = -Kn-l{z) - -Kn{z) = -Kn+l(z) + -Kn{z), 

k'„{z) = -zKa^iiz), k'^{z) = -Zi^„+i(z). (V.IO) 
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^^11 1u'j^-um db , {u' ■b-u-b'y)u'j^ + {u-b-u' ■bj)uM\ nm\ 

~-T\'^M + /„. ,./xo„.o . (V.ll) 



Using 

dW ~ TV ' ~ b V"'" ' {u ■ u')^v^ 

one is led to the final result 

(•27r)"+i ( bz - - \ 

= - UM)Kr,{z) + ibMKn+l{z)j • (V.12) 

(c) The integral 

J {q ■ u)^ J [qouo - qwyv)^ q^ 

after integration over dq^5{q^) gives the extra factor (— q • u)""^, which by virtue of rf(q • u)6{ku + q • u)/'yv 
gives l/{k ■ u)^. Thus 

Similarly, for the vectorial integrals 

riA)_ [ n Aq-u')5{{k-q)-u)e^'i-\ _ Im ^ 
=jd Q ^^7^2 IM - (V.15) 

Similar relations hold for the primed integrals. The denominators are (A; • u')"^. 

(d) In the case of Adi x T'^ space-time, one is led to the invariant integral 

f S{q-u')6{{k-q)-u)e^'^->> 



where n = vn^ and n = 27rn/L, n e Z . It is computed in the same way. After tlie g -integration, the 
denominator becomes q^ + li^ + (fc • u/'fv)'^, and leads to 

In = Ko{Z^), (V.17) 



with Zn = V z'^ + n^6^ and z = k ■ u b/{'yv). If m represents the single one-dimensional KK mass mi, we will 
refer to Zi = \J -\- mfb"^ and to the integral 7; respectively. 

(e) The vectorial integral can be computed by differentiation and the result is 

Analogously, if the integrals have the factor in the denominator {q ■ uY of the integrand, give the additional 
factor (fc • u)~^ by virtue of the delta-function 6{{k — q) ■ u), namely one obtains 

=I^(k. = I'M (fc . u)-\ (V.19) 

Analogous relations hold for the primed integrals. 

VI. APPENDIX II 
A. Angular integrals 

(a) In the main text the following angular integrals 6 were needed for integer m and n 

7r 

VZ = J^-^d0, ^ = l-vcose. (VI.l) 
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Making use of the formula [25|, valid for any real a > and Rei^ > 0: 



de^{~-] ^T{.){a^ - b^r-^y^-^'^Pli'-^,, [ ) , (VL2 







we express the result in terms of the associated Legendre function Pj^{z). In our case a = 1, = — w, so: 

I/- 1/2 / „\ n/2 



y2u-l 



(') " ^n-)l-^'^'''Pl'-;-\,2 (7) = (')' ^/^^ (!^) T-Z^+^/^^/f ,(7). (VI.3) 



For 7^1, one can use the asymptotic formula |24j : 
For 2m > n + 1 one finds to leading order 
while for 2to < n + 1 

on-m p / n+1 A p z' »+l _ ^A 

" W2 r(n-m+l) ^ ^ 

In the case 2m = n + 1 an expansion of the integral is logarithmic, 
(b) Another integral used in the main text is: 

|sinV.rd»= ^, z._2+iv A' ■ (VI-7) 

^ ^ \ 2 I 

B. Integrals of products of Macdonald functions 

Computation of integrals over the frequency or the impact parameter involving products of two Macdonald 
functions of the same argument is performed using the formula [25| : 

KJcz)K,{cz)z"-'dz = ^ 2 >^ ^ \^ ' ^ ^ ' A (VI.8) 

u 

Actually, only such integrals are needed in this work. More general integrals involving functions of different 
arguments arise from interference terms. They can be computed using the formula 

CO J. 

i^.(Az) = A''5]-^(l-A2)'=i^,+fc(z), |1-A2|<1. (VI.9) 

fe=0 

The typical integral 

oo 

e = y (VI.IO) 



where z — can be cast for 7^ < 1 (i.e. ^ 6* < ^2/7) into the form 

(70"(y) j K,{i^z')K,,{z')z'^+^'z'^-^dz', (VI.ll) 



and for 7^ > 1 {^JVl ^ ^' =^ tt) to 



1 /7W 



(7O™' V P 



D-l 



K^{z)K^, (z/70z'"+"''z-°-2dz. (VI.12) 
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One is interested in an estimate of (jVI.10[) in powers of 7 ^ 1. Use (|VL9p to write 

D-l 00 

-] 

P 



00 



z'^'-^dz', 7^ < 1; 



fe=0 
D-1 00 



(70 



^[rpj J ^'K..{z)K.,^,iz)z-'^-'z^-'d, 7e>l. 



(VI.13) 



Denote by M = m + m! + — 1 and integrate over z with the help of fVLS)) . The integral in (|VI.131 a) is equal 
to 



(VI.14) 



while the one in (j VI. 131 b) is related to the above by the exchange v -ir^ v' . 

Given that for fixed A > the function r(x)r(A — x) is decreasing for x in [0, A/2] and increasing in [A/2, A], 
one obtains the inequality 



(VI.15) 



which leads to the estimate 



e < 



pz-sp (M+p/^ p (m^^ (^^)-(™'+.') [v/£^py-- f (1 - ^) r(:.' + k)/k\, 7C 
Both series are summable because 

r(n) 



< 1; 
> 1. 



$:l(i-aYr(n + fc) = ^, H<1, 



fc! 

fe=0 

and give the following upper bounds, which are enough for our purposes: 

D-l 



(VI.16) 



e< 




(M±^)p( 



M-v' -V 



7^> 1. 



(VI.17) 



For V = v' the estimate is continuous at 6* = arccos((7 — l)/7'i;). 



VII. APPENDIX III 



A. Effective number of interacting $ KK modes 



In this subsection we would like to count the effective number of contributing interaction KK modes of the 
bulk field $, studied in Subsection III.B, and estimate the error due to the replacement of the summation over 
those by integration. 

Start with the case D — 5, i.e. with one-dimensional mass sequence. It was argued that the dominant 
contribution to the energy loss comes from the first term in the parenthesis of (|III.20|) and requires the evaluation 
of the two-parameter integrals 



G(a,/3)= / Ko(Vz^ + a^) Ko(Vz^ + P^) z^dz, 



(VII.l) 



for a ^ 0, /? ^ and a + (3 > 0. Using the integral representation [24 



2 



K^{x)K^{y) = - / exp 



2t 



(VII.2) 



we obtain the following approximation G(a, /3) for G: 



(VII.3) 
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In view of the exponential fall-off of G it is natural to separate the modes to light for a + (3 < 2 and heavy for 
a + /3 > 2. The choice a + (3 — 2 for the boundary is justified by just looking at the numerical plots. For the 
heavy modes one can neglect the unity in the parenthesis of (|VII.3P and write 

G(a,/3) « ^-'^'e-i^^'^J^^^ , " + > 2 (VII.4) 
Correspondingly, for the light modes one may take approximately 

G(a,/3)«^e-("+^), a + /3 < 2. (VII.5) 

In the cases of interest here a = 27Tbl/L,(3 — 2Trbl'/L with integer I and I'. Thus, denoting hy K = I + I' , one 
sees that the integer part 

= [L/{TTb)] (VII.6) 

defines the boundary between light {K < Kc) and heavy (K > Kc) modes. If L < irb, there are no light modes 
in this classification. In the sequel only the case L > nb {Kc > 1) will be discussed. 

To estimate separately the contributions of the light (heavy) interaction modes to the emitted energy, sub- 
stitute the first term of pil.20p into (III.16P and integrate over angles, using (IVII.5|) and (|VII.4p . respectively. 
For the light modes one obtains (up to coefficients) 

E flight ^^T. e-'^^'/^(i^ + 1), Kc ^ 1. (VII.7) 
It is easy to sum the geometric series for any Kc- In particular, for Kc ^ 1 the highest power in Kc gives 

Y^Gu^^^c^^^Kl Kc»l, (VII.8) 
with Alight = 1/4 - 3/(40^) « 0.148 « 1/7. Similarly, the heavy modes give 

_3/2 -2K/K, K _3/2 °^ 

K=K^ + 1 l=Q l^^^^c K=K^ + 1 

Replacing the sum over K by integration from Kc to oo, one obtains for large Kc'. 

E Ghoavy - ^^f^Kl Kc » 1, (VII.IO) 

with Aheavy = 7 / (v^e^) + 30F [l - $ (%/2)] /4 w 0.73 w 5/7, where $ is the Laplace's error function. In the 
case L < nb (no light modes) one should integrate from to cx) and substitute Kc — t- L/{TTb) in (|VII.10p . 
Using (|VII.8[) and (jVII.10[) one obtains the following ratio for d = 1: 



A-^heavy 2"'^'^^Aheavy 

Thus, despite the fact that each heavy mode is exponentially suppressed their total contribution is comparable 
to the one of the light modes. 

The above generalizes to arbitrary d > 1. One has to integrate over the domain W defined as Af = fi 4- fi' < 
2ttKc/L = Mc which represents the multidimensional bispherical octahedron (the coefficient 7r^/32 is omitted): 

E4?M - ^^1^ / e-^^/^^ifrirr-^cindii' ^ ^B(d, d)7(2d, 2) Xf, (VII.12) 

w 

with j(x,y) is the lower incomplete gamma-function, B(a;,?/) is the Euler beta-function. The corresponding 
heavy mode contribution reads: 

^^(.) 2V^(yA_i)^ /- e-^^/^- ,__,,,^_^_, »LiB(d+l,d+l)r(2rf+l/2,2) 

w 

with r(x, y) the upper incomplete gamma-function. Thus, the ratio is 

SGght ^V^B(d,d)7(2d,2) 

y-gW - 2V2B(d+l,d +l)r(2d+ 1/2,2) ' ^ ■ ' 

^ heavy 
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a decreasing function of d. 

Add (|VII.12|) with (jVII.lSp for the emitted energy in _D— dimensions and compare to the energy loss AE4 in 
the purely four-dimensional case to obtain the estimate 



AE4 



Kf. (VII.15) 



This is in line with Table I and explains the origin of the enhancement in the entries of the first row (bulk 
interaction dominance), as compared to the ones in the second row (the case of brane interaction dominance). 



B. Effective number of emission modes and the angular distribution 

Let us compare the role of KK emission modes, studied in Subsection III.C, with that of KK interaction 
modes analysed in the previous above. Passing to integration over the emission modes, we have found that for 
ultrarelativistic velocities the angular-frequency distribution is the same as in the non-compactified case: the 
typical frequency being (uj) = 7^/& and the emission angles within the cone < = arccosu. From equation 
(IIIL34P one can see that only a finite number of emission KK modes contributes to the total energy loss. Thus, 
unlike the case of interaction modes, here one can determine the effective number of emission modes from the 
energy loss right from the beginning. Indeed, from v < cos a < 1 and sin a < 7"-*^, one concludes that 

n < wj^^ = 27/6, 

so that the effective number of emission modes is 

^cff = (VII.16) 

This is 7 times greater than the effective number of interaction modes iVoff = jKc- Notice that the massive 
arguments (with equal masses li of the scalar quantum) of the Macdonald functions in the rest frame of the fast 
particle, do not depend on the emission angle: 

Zn = -f n2i;2, A^^/2 + n2, (VII.17) 

V 

with uj' = |k| in this frame. From (jVII.17[) it is clear that the effective number of modes (demanding the 
argument to be not greater than unity) of z'™) is times larger than the corresponding one of Z^. The overall 
effect of massive emission modes is found passing to integration assuming that N^e is large. 

The angular distribution for a specific KK-mode is more complicated since z*^^^ depends on uj,N and the 
angle 0. For light emission modes one has cos a > v. The form ofE — l~v cos a cos suggests the introduction 
of the "effective velocity" 



V COS a, 



which differs from unity by 0(7 ). Thus, the angular and spectral properties can be derived from the effective 



quantities w^g ' and 7^^' 



-1/2 



AN) 



For the light modes one can expand 



(VII.18) 



274' ^^-""-s 2vj^' 

^cff -^-TT— 4' %ff -7 1-— ^ • (VII.19) 



Thus, the correction is of the order of 0(7 2) and one can substitute the frequency by its average value: 

?eff COS^y^g -COS& 

2^74' '°ff ' 72 

This leads to the effective replacement of the emission angle 

0is^ = yj0^ + N^ b^/v-f^. 
Therefore the emission of a given N-mode is concentrated within the cone 



= ^0"^ + b'^/v-i^ = arccosw^ff \ (VII.20) 
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and the small-angle approximation is valid for all light modes. 

The radiation flux of a single KK-mode N under massless interaction mode can be obtained substituting the 

effective Lorentz factor 7^^'* into Eqn. (jll.46[) (with D = 4). Note that although it is two powers of gamma 
smaller than the general result, the exact correction can not be calculated by such analysis, because of the 
presence of the interference term with the product of two Macdonald functions with different arguments. The 
exact results may be obtained in the center-mass frame, where the interference term vanishes. 

Total radiation flux can be estimated simply as the product of the four-dimensional result by the total number 
massive emission modes: 

E4+d — EiN^^ . 

Again, this explains the origin of the relative enhancement in the entries of the first column (corresponding to 
bulk emission) of Table I, compared to the ones of the second (radiation of the brane field). 
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